INTEGRABLE GL(2) GEOMETRY AND HYDRODYNAMIC PARTIAL 

DIFFERENTIAL EQUATIONS 

ABRAHAM D. SMITH 

Abstract. This article is a local analysis of integrable GI/(2)-structures of degree 4. A GL{2)- 
structure of degree n corresponds to a distribution of rational normal cones over a manifold of 
dimension 71 + 1. Integrability corresponds to the existence of many submanifolds that are spanned 
by lines in the cones. These G'L(2)-structures are important because they naturally arise from a 
certain family of second-order hyperbolic PDEs in three variables that are integrable via hydrody- 
namic reduction. Familiar examples include the wave equation, the first flow of the dKP hierarchy, 
and the Boyer-Finley equation. 

The main results are a structure theorem for integrable GL(2)-structures, a classification for con- 
nected integrable GL(2)-structures, and an equivalence between local integrable G-L(2)-structures 
and Hessian hydrodynamic hyperbolic PDEs in three variables. 
00 This yields natural geometric characterizations of the wave equation, the first flow of the dKP 

hierarchy, and several others. It also provides an intrinsic, coordinate-free infrastructure to describe 
a large class of hydrodynamic integrable systems in three variables. 
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Introduction 

A fundamental problem in analysis is to understand why some differential equations (particularly 
hyperbolic equations arising from wave-like equations or from differential geometry) are integrable. 
This problem is compounded by the many competing notions of integrability [35]. Informally, 
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a PDE is called integrable if it can be solved thanks to the existence of an infinite hierarchy of 
conservation laws or of a family of invariant foliations along characteristics. Integrable PDEs of 
all types are well-studied in two (1 + 1) independent variables; however, examples are increasingly 
rare and poorly understood in higher dimensions. 

A particularly interesting class of integrable PDEs are those that can be integrated by decom- 
posing the equation to a set of coupled first-order equations that represent traveling waves [3T] . 
This technique is called "integration by hydrodynamic reduction," and it has been extensively 
studied for various special classes of second-order PDEs in three or more independent variables 

El [la El [321 E3] 

Consider a second-order hyperbolic equation of the form F{^\ u, Ui, Uij) = on scalar functions 
In the method of hydrodynamic reduction, one hopes to integrate the PDE F = 
by stipulating that the solution function u (and its derivatives) may be written as u{R^, . . . ,R^) 
for an a priori unknown number k of functions . . . , whose derivatives admit the 

"commuting" relations 

d ■ d ■ d ■ d ■ 

(1) -g^R' = pliR)-^R\ -g^R' = pI{R)^R' 

which also imply 

(2) 1 d^^ 1 

^' Pl-pidR^ Pl-pidR^^ 

Equations ([T]) and ^ can be solved by other methods [31]. The reduction of F = to Equation ^ 
is called a fc-parameter hydrodynamic reduction. The equation F = is called integrable via k- 
parameter hydrodynamic reduction if it admits an infinite family of fc-parameter hydrodynamic 
reductions and this family is itself parametrized by k functions of one variable. The functions 
are called Riemann invariants, and their level sets define the foliations along char- 
acteristics in the informal definition above. See the works cited above, particularly [H], [15] and 
[T7] . for examples and detailed exposition of this technique. 

Many examples of these so-called hydrodynamic equations are known, such as the wave equation 
U22 = Ui3, the dispersionless Kadomtsev-Petviashvih (dKP) equation U22 = {ui — ^^3)3, the first 
flow of the dKP hierarchy U22 = M13 — 1(^33)^, the Boyer-Finley equation U22 + M33 = e"", and 
some well-known integrable hierarchies. Tests for hydrodynamic integrability also exist, but so far 
there is no intrinsic, coordinate-free theory that classifies and generates all of these equations. 

For the case of second-order hydrodynamic equations in three variables that involve only the 
Hessian of the solution function, recent work by Ferapontov, et ai, reveals a relationship between 
the integrability of the PDEs, the symplectic contact symmetries of the equations, and the equa- 
tions' underlying degree-4 GL(2)-structures [H]. Integrable G'L(2)-structures are not well-known 
in the literature; they arise in the search for exotic affine holonomies in differential geometry [5] 
and a similar conformal geometry appears useful in the equivalence problem for ODEs [24], but 
the application to PDEs is very recent [I^ E|. 

The purpose of this article is to fully analyze the local geometry of integrable GL(2)-structures of 
degree 4, with an eye towards understanding PDEs that are integrable via hydrodynamic reduction. 
The main result is a natural classification of local integrable G'L(2)-structures and the associated 
PDEs. This is achieved using geometry d la Cartan, namely exterior differential systems, moving 
frames, and the method of equivalence. 

Here is a summary of the contents of this article, including abbreviated and non-technical ver- 
sions of the main theorems: 

Section [1] defines G'L(2)-structures and their 2-integrability and 3-integrability. The property of 
fc-integrability involves the existence of "many" foliations by certain submanifolds of dimension k. 
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Section [T] also introduces the motivating observation by Ferapontov et al. that hyperbohc PDEs in 
three variables that are integrable by means of hydrodynamic reductions lead to G'L(2)-structures 
that are both 2-integrable and 3-integrable (abbreviated 2,3-integrable). 

Section [2] contains a description of the necessary SL(2) representation theory and an explicit 
solution to the first-order equivalence problem for GL(2)-structures. Corollary 2.2 defines a global 
preferred connection for GL(2)-structures that has essential torsion in the direct sum of M.^, M7 , 
MP, and M^^, each of which is an irreducible SL{2) module. 

In Section |3| the exterior differential system describing 2-integrability for a GL(2)-structure is 
analyzed. The most important object in this article, the irreducible M^-valued torsion, is first 



emphasized in Theorem 3.2 



Theorem |3.2[ A GL (2) -structure is 2-integrable if and only if its torsion takes values only in 



Section |4| provides Theorem 4.3, a complete local description of GL(2)-structures that are 2,3- 
integrable. The immediate consequence of this theorem is that the value of the torsion at a single 
point completely determines the local structure of a 2,3-integrable GL(2)-structure. 



Theorem 4.3, A GL (2) -structure is 2,3-integrable if and only if its structure equations are of the 
form 

du = -9 Auj + T{uj A u) 
(3) <\e = -6 AO + T'^{uj Au) 

dT = J{T){uj,e) 

for an essential torsion T in M.^ and a 9 x 9 matrix J{T) depending on T , as in Appendix p4| 

Section [5] presents a classification of connected 2,3-integrable G'L(2)-structures. This classifica- 
tion is provided by a natural stratification of into equivalence classes that are defined by the 



matrix J{T). These equivalence classes are explicitly identified in Theorem 5.3 



Theorem 5.3[ Every connected 2,3-integrable GL(2)- structure belongs to exactly one of 55 equiv- 
alence classes, which are given by a stratification of into the root factorization types of real 
binary octic polynomials. 

In Section |6l the structure equations from Section |4] are applied to reproduce the integrable 



PDEs, yielding Theorem 6.3 



Theorem 6.3, Every local 2,3-integrable GL{2)- structure arises from a Hessian hydrodynamic 



hyperbolic PDE 



Therefore, the classification in Theorem |5.3 is also a classification of Hessian hydrodynamic 
hyperbolic PDEs. Some previously known examples of such PDEs are identified in the new clas- 
sification, and several new integrable PDEs are constructed. Also, a relationship between Hessian 
hydrodynamic hyperbolic PDEs in three variables and hyperbolic PDEs in two variables (Monge- 
Ampere, Goursat, or generic type) is observed. 

This work would have been impossible to complete without the uncanny speed and accuracy 
of computer algebra software, and there is no sense in retyping dozens of pages of formulas and 
identities. Thus, a reader looking for explicit computational details may be disappointed by this 
article (in which case I direct the reader to [27]). However, most proofs rely only on basic concepts 
from differential geometry, linear algebra, and representation theory, and a reader comfortable with 
the standard methods of exterior differential systems (and their various software implementations) 
can readily supply the computational details. A Maple file containing the structure equations 



appearing in Theorem 4.3 and the Maurer-Cartan form appearing in Theorem 6.3 will remain 
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available at flE] as long as possible. This file should also be bundled with this article on |ar Xiv.org 



The content of this file allows rapid reproduction of all results following Theorem 4.3 

Please note that the terms "hyperbolic" and "integrable" are heavily used in slightly different 
contexts throughout this article. "Hyperbolic" refers to hyperplanes having maximal intersection 



with certain projective varieties (for example, in Lemma 1.2), to PDEs in both two and three 
independent variables with appropriate leading symbol (for example, the wave equation), and to 
exterior differential systems with appropriate tableaux (for example, in the proof of Theorem |3.1[ ). 
The term "integrable" refers both to PDEs that admit infinitely many exact solutions and to GL{2)- 
structures that admit many secant submanifolds (submanifolds which one might be tempted to 
also call "hyperbolic"). Of course, all of these concepts are intimately related, so the standard 
overloading of these definitions is reasonable. 

I owe many thanks to Robert Bryant and Niky Kamran for their wisdom, guidance, and support, 
to Dennis The and Eugene Ferapontov for their many helpful comments, and to Jeanne Clelland, 
whose "Cartan" package for Maple made the computations bearable. 

1. Background 

Let Vn denote the vector space of degree n homogeneous polynomials in x and y with real 
coefficients. Identify V„ with M""*"^ using the terms from the binomial theorem to produce a basis; 
for example, V4 — by 

(4) f_4 x'^ + f _2 ix^y + Vq 6x^2/^ + f 2 4:xy^ + f 4 y'' ^-7■ (f_4, f _2, Vq, V2, v^) G M^. 

Let M denote a 5-dimensional smooth manifold, and let J-" denote the V4-valued coframe bundle 
over M, so the fiber J^p is the set of all isomorphisms TpM — )■ V4. The coframe bundle J-" is a 
principal right GL{V4) bundle. 

Recall that V„ is the unique irreducible representation of st(2) of dimension n+1. The action of 
s((2) on Vn is generated by 

d d d 

(5) X = y—, Y = -X—, and H = x- y—. 

ox oy ox Oy 

Definition 1.1. A GL {2) -structure B ^ M is a reduction of T with fiber group GL(2) C 
GLiVn) infinitesimally generated by X , Y , H , and the scaling action I . 

Note that throughout this article, the symbols GL{2) and 0t(2) are used to denote the abstract 
Lie group and its Lie algebra as well as any particular irreducible representations GL{2) C GL(Vn) 
and 0t(2) C 0t(V„). Since SL{2) has a unique irreducible representation Vn in each dimension n+1, 
it is always clear from context which representation of GL{2) is being considered. 

Generally, a GL(2)-structure is said to have "degree n" if the base space M has dimension n+1. 
The degree 3 case is thoroughly studied in [5], and various observations are made for all n in [27] . 
This article considers only the degree 4 case. The PGL{2) action generated by X, Y, and H is 
the symmetry group of a rational normal curve in P(V4). A rational normal curve is a curve of 
degree n in P(]R"'"''^). All such curves are PGL{n + l)-equivalent to {[g"^ : g"'~^h gh^'^ : (7"]} 

j20] . The de-projectivization of the rational normal curve is the rational normal cone, which has 
a GL{2) symmetry group and is usually described as {{g'^, g'^~^h, . . . ,gh^~^,h"')} C IR""'"^ Thus, 



the geometric content of Definition 1.1 is contained in the following lemma. 



Lemma 1.1. A GL (2) -structure B — )■ is equivalent to a distribution of rational normal cones 
C C TM. For any b e Bp, b{Cp) = {{gx + hyf : g,h eM.} <Z V4. 

Definition 1.2 (Integrability). Let B M be a GL{2)- structure. A k- dimensional linear sub space 
E C TpM is k-secant if EnCp is k distinct lines. Equivalently, E is k-secant if E is spanned by 
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k vectors in Cp. A k- dimensional suhmanifold N C M is k-secant ifTpN is a k-secant subspace 
ofTpM for every p & N. A GL{2) -structure is k-integrable if, for every k-secant linear subspace 
E G Gri^{T!M), there exists a k-secant submanifold N with E = TpN. 

Locally, 1-integrability is uninteresting, as it is simply the existence of a local flow for a vector 
field. Since M has dimension 5, being 4-secant is an open condition on G Gr^lTM); however, 
the condition is closed for k=2 (bi-secant) and k=3 (tri-secant). A GL(2)-structure that is both 
2-integrable and 3-integrable is called "2,3-integrable," or simply "integrable." 

Notice that all the definitions are projectively invariant. Indeed, one could equivalently study 
the principal right PGL{2) bundle over M defined by the symmetries of rational normal curves 
in PTM. When considering the de-projectivized GL[2) geometry, one occasionally encounters an 
additional Z/2Z symmetry by ±/. For the purposes of this article, it is easier to proceed using 
affine geometry and deal with this symmetry when it arises. 

The geometric content of this article is a complete local description and classification of GL{2)- 
structures of degree 4 that are 2,3-integrable. Such G'L(2)-structures are particularly interesting 
because of their deep connection to integrable PDEs. The remainder of this section outlines how 
Ferapontov, Hadjikos, and Khusnutdinova obtain 2,3-integrable GL(2)-structures from certain 
hydrodynamic PDEs in |14j . 

1.1. Hessian Hydrodynamic Hyperbolic PDEs. The space of 1-jets over — t- M is = 
{(^*, z,pi)} where I < i, j < 3. The space of 2-jets is = {(^-', z,pi, Uij)} where 1 < i,j < 3 and 
Uij = Uji. The contact system on is the differential ideal J'^ generated by 

(6) {dz-pAC, dp,-U,,de}. 

A second-order PDE on M : ]R3 _^ M is equivalent to the level set E ^(0) of some smooth function 
F : — )■ M on which the projection — )■ is a submersion, and a solution is a function 
u : M'^ — )► M whose jet-graph is a subset of E~^{0). 

Recall that a section of is the jet-graph of a function if and only if it is an integral of the 
contact system. A contact transformation on is a diffeomorphism : — ?■ such that 
ip*{J'^) = and one studies PDEs geometrically by examining the properties of -F^^(O) that are 
invariant under contact transformations. 

By a theorem of Backlund Theorem 4.32], every contact transformation on is the prolon- 
gation of a contact transformation on J^. Therefore, for any contact transformation there exist 
functions Aij, Bij, Cij, Dij, ai, rrii, ai, rhi, and c, (written as four matrices, four column vectors, 
and a scalar) such that 

/dA fBmC\ /dA 

(7) 7/-* d2 = m* c a* d^ . 

\dpl \A a dJ \dpj 

The assumption that is a. contact transformation implies that the symplectic form a = dpi A 
d^^ -|- dp2 A d^^ -|- dp3 A d,^'^ is preserved up to scale, so ip*{(y) = Aa, A 7^ 0. One can compute 
Uij = 'il)*{Uij) by noting that il)*{(\.pi — Uijd^^) = modulo J''^, so in matrix notation 

^p*idp) - t/^*(dO = (Ad^ + adz + Ddp) - U{Bd^ + rhdz + Cdp) 

(8) 

= {A + ap' + DU)d^ - U{B + mp* + CU)di. 
In particular, under the non-degeneracy assumption ip*{dS} A d,^^ A d,^'^) 7^ 0, this yields 

(9) tj = {A + ap' + DU){B + mp' + CU)-\ 
For this article, consider only PDEs on -u : M'^ — )■ M of the form 

(10) E{uii, U12, Ui3, U22, U23, M33) = 0. 
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Instead of considering the orbit of this PDE under all contact transformations, consider only those 
contact transformations that preserve the class of Hessian-only equations, those of the form in 



Equation (10). That is, consider only those contact transformations ip such that ilj*{dUij) = 
modulo {dUki}- Then it must be that a = m = in Equation ([t]). Thus, the contact trans- 
formations that preserve the Hessian-only form of Equation ( [Io| are elements of the conformal 
symplectic group, 

CSp{3) = {g eGL{6,R) : a{gv, gw) = \a{v,w) yv,w {any \ E R*)} 

(11) f fB C 
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= A'B- B'A = D'C - C'D, Xh = D'B - C'A 



In the most general case, the completed domain of all Hessian-only F's is the Lagrangian Grass- 
mannian, 

(12) A = f/f = {UeGr3(R^) : (t\u = 0}, 
with local coordinates given, for example, by its non-,^-vertical open subset 

(13) A° = {f/G Aid^^Ad^^Ad^^^lf/T^O}. 

The conformal symplectic group forms a bundle, H : CSp{3) — )■ A. In fact, when considering 
local transformations of PDEs, only the transformations in CS'p(3)° = H~^(A°) are permissible, as 
the condition d^^ A d^^ A d^^ 7^ must be preserved. As above, one could eliminate the conformal 
factor and consider only 5*^(3)° actions by projectivizing, but there is little utility in doing so. 

To be explicit regarding the bundle, suppose U G A°. The open manifold A° has coordinates 
given by the components of symmetric matrices, Uij = Uji, and dpi = Uijd^^ modulo J'^. Then 

(") *' (^0 - {1 i) (^^ 

The assumption that d.^ and ?/'*(d^) both have maximum rank when pulled back to implies that 



(B + CU) is nonsingular. Moreover, Equation (14) shows that CSp{3)° acts on A° by 

(15) g : U ^ g{U) = {A + DU){B + CU)-\ 
The fiber subgroup over U G A° is the stabilizer of U, so 

(16) U-\U) = {g:U = g{U) = {A + DU){B + CUy'} ^{g:A = ^}. 

Because Jg is in the fiber over G A", the projection H : CSp{'i)° A° can be computed from 



Equation (15) as 

(17) Ii:g^ (7(0) = AB'K 



The most important fact about Equation (10) is that M = F ^(0) admits a natural GL{2)- 



structure wherever the PDE is hyperbolic. Recall that a second-order PDE is called hyperbolic if its 
leading symbol is a non- degenerate matrix with split signature, and this property is invariant under 
contact transformations. A related notion of hyperbolicity for hyperplanes sitting in projective 
space is needed to describe the induced G'L(2)-structure. The complex Veronese variety is the 
2-dimensional projective variety 

(18) {[ZiZi : : : Z^Z^ : Z^Z^^Z G CP^} C CP^ 

The intersection of a generic hyperplane with the Veronese variety is a 1-dimensional rational 
normal curve, and this curve is uniquely given by such an intersection. The real case needed 
for the present PDE theory requires a bit more detail to describe accurately. Consider MP^ with 
coordinates [W^i : ■ ■ ■ : Wq\^ MP^ with coordinates [Zi : Z2 ■ Z^], and the Veronese variety defined 
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over M as above. A generic hyperplane in MP^ is defined uniquely (up to scale and sign) by a single 
equation 

(19) anWi + ai2W2 + aisW^ + 022^^4 + 023^^5 + 033^6 = 0. 

The intersection of this hyperplane with the Veronese variety yields an equation on M^: 

(20) aii(Zi)^ + 012^1^2 + ai3^i^3 + ^22(^2)^ + ^23^2^3 + ^33(^3)^ = 0. 



Depending on the (real) coefficients aij, Equation (20) may or may not have real solutions Z. If 
Equation (20) has real solutions, then the de-pro jectivized solution in is a real quadric surface. 



This quadric may or may not be degenerate. The existence of real nondegenerate solutions is an 
open condition on the hyperplane in the topology of (j'r5(M^). If this condition is satisfied, the 
hyperplane defined by {atj} is called hyperbolic. The precise algebraic condition for hyperbolicity 
is that the real symmetric matrix (aij) is nonsingular and has split signature Using the 

identifications = Gr^iR^) = Sjm^{R^), the cone over the Veronese variety (Veronese cone) 
is identified with the space of symmetric matrices (a^) such that rank(ajj) < 1. A hyperbolic 
hyperplane in intersects the Veronese cone in a rational normal cone, and every rational normal 
cone in can be written (uniquely, up to scale) this way |20j . 

Lemma 1.2. Let F : A° — > M 6e a smooth function, and suppose U E A° such that F{U) = 0, 
dFu 7^ 0, and ker(dF[/) is hyperbolic as a hyperplane in TjjA" = Sym^(M^). Then there is an 
open 5-dimensional submanifold M C A° defined by F\m = in a neighborhood of U, and M 
admits a distribution C of rational normal cones. That is, M admits a GL{2)- structure. 

Proof. Aside from the condition on hyperbolicity of the tangent space, the lemma is simply a 
statement of the implicit function theorem. Here is an explanation of hyperbolicity and its relation 
to G'L(2)-structures. 

Fix a coframe s on A° such that su '■ T^yA" — i- Sym^(]R"). (An obvious choice is s = d?7, using 
the symmetric matrices as coordinates for A°.) So, there is a distribution of 3-dimensional Veronese 
cones defined by {P G TA° : ranks(P) < 1}. This distribution of cones is CS'p(3)-invariant, for 
in the first part of [9l Theorem XX], Cartan proves that the above bundle CSp{?)) is exactly the 
frame bundle over A whose action on T^/A ~ Sym^(M") is to act irreducibly on the coefficients of 



the equation of a cone, for example Equation (20). 



In the context of the lemma. Equation (19) describes the intersection of ker(dF[/) with the 
Veronese cone, so the symmetric matrix {aij) is, up to scale, the leading symbol of F at U . 
Therefore, if [/ is a regular point for the regular value of F : A° — M and if the PDE F = 
is hyperbolic at U, then F~^(0) admits a field of rational normal cones near U given by these 
intersections. □ 



The fact that every hyperbolic PDE of the form in Equation (10) admits a GL(2)-structure 
would only be an algebraic curiosity, except that the integrability of the PDE is intimately related 
to the 2,3-integrability of the GL(2)-structure. 



Theorem 1.3 (Theorem 3 in [H]). Fix a hyperbolic PDE of the form in Equation (10) and 
its corresponding GL (2) -structure tt : S — )■ M. Then B is 2-integrable. Moreover, the PDE is 
integrable via 3-parameter hydrodynamic reductions if and only if B is also 3-integrable. 



Theorem 1.3 as presented in [T4j apparently requires that, for each 3-secant A^^ C M^, the char- 



acteristic net defined by the intersection TA^ fl C is a coordinate net. According to Corollary 4^ 
this requirement is redundant. 

As discussed in [T3], for PDEs in three independent variables, the existence of 3-parameter 
hydrodynamic reductions parametrized by three functions of one variable implies the existence of 
/c-parameter hydrodynamic reductions parametrized by k functions of one variable for all A; > 3. 
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Corollary 1.4. F{uij) = is integrable via k-parameter hydrodynamic reductions for all k > 2 if 
and only if the induced GL (2) -structure over F^^{0) is 2, 3 -integrable. 

Definition 1.3 (Hessian hydrodynamic PDEs). A PDF on u : M.^ is Hessian hydrody- 

namic if and only if it is hyperbolic, is of the form F{uii, . . . , M33) = 0, and is integrable by means 
of 3-parameter hydrodynamic reductions. 

Tlie remainder of tliis article is a study of Hessian hydrodynamic PDEs via their associated 
2,3-integrable G'L(2)-structures. 

2. A Preferred Connection 

In this section, Cartan's method of equivalence is applied to GL(2)-structures. Cartan's method 
of equivalence is a standard tool in the field of exterior differential systems; it is an algorithm 
for finding all of the differential invariants of a geometric structure, and the first step is to fix a 
preferred connection among all the possible pseudo-connections of a geometric structure [U |19|, 
[22] • (The distinction between a connection and a pseudo-connection is simply whether it was 
obtained via such an algorithm.) The result of the algorithm is a preferred global coframe for the 
GL(2)-structure that splits into the semi-basic "tautological" form, u, and the vertical 0[(2)-valued 
"connection" form, 6. 

For a GL(2)-structure vr : i3 — )■ M, let u denote the tautological 1-form defined by Ub = b o n^, : 
TbB — i- V4. As always, a; is a globally-defined 1-form on B with linearly independent components, 
and it is semi-basic, meaning wlkcrTr, = 0. A pseudo-connection on i3 is a 1-form 6 taking values in 
the non-trivial representation of st(2) in 0[(V4) such that du = —9/\uj + T{uj/\uj) for some torsion 
T{h) : V4 A V4 ^ V4. 

The goal is to fix a particular 9 that minimizes the torsion T. To carry out this process, one 
needs to understand the SL{2) representation theory on V„. The decomposition of the tensor 
product Vm ® V„ into irreducible components is 

(21) Vm®Vn = V 

The projections onto the various components are given by the Clebsch-Gordan [51 [21] pairings 
(■, ■) : Vm © V„ — Vm+n-2p, which are provided by the formula 



(22) Mp = ^^h 



p! \k) dxP-^dy^ dx^dyP-^' 

This pairing has some important properties. Notice that {u,v)p = {v,u)p and that the 

pairing is nontrivial for < p < min{m,?2}. Hence, the tensor decomposition can be further 
refined in terms of the symmetric and alternating tensors: 

V„ O V„ = Vsn © V2„_4 © ■ ■ ■ © Vo or 2, 
^ ^ V„ A Vn = V2n-2 © V2„-6 © " " " © V2 or 0- 

Notice too that (■, ■)^ : Vn © Vn — Vo = M is a non-degenerate symmetric or skew bilinear form. 
Hence, for fixed u EVn the map {u, ■)^ : V^ — )■ Vq = provides a natural identification, Vn = V*, 
and one need never distinguish between dual spaces when considering representations. 

For any derivation over ]R[x, y], a Leibniz rule over the pairing holds. Because SL{2) is infinites- 
imally generated by X, Y, and H, this means that the pairings are 5'L(2)-equivariant. That is, 
a{{u, v)p) = {a{u),v)p + {u, a(f ))p for any a G s((2), which implies that a ■ {u, v)^ = {a ■ u,a ■ v)^ 
for any a G SL{2). The pairing is not GL(2)-equivariant, but the scaling action is easily computed 
easily where required. The geometric objects encountered here are projectively defined, so this 
variance in scaling is of little concern. 
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Most importantly, the pairing can be generalized to binary-polynomial- valued differential forms 
over a manifold. If m G r(A''T*B V^) and v e r(A^T*B Vn), then extend the definition by 
using the wedge product: 



(24) 



{u 



A 



QxP-^dy^ dx^dyP~^ 



In this gencrahzation, the symmetry of the pairing is further altered by the degree of the forms: 
{u, v)p — (— 1)^*+P (i), u)p. If A is an R-valued 1-form acting by the scaling action Xl^n on Vm, then 
X Au may be written as the trivial pairing (A, u)^ = (— 1)^' (m, A)q G r(A''"'"^T*,B (g) Vm)- 

As SL{2) representations, 0t(2) = si{2) © M = V2 © Vq. Thus, a 0[(2)-valued pseudo-connection 
9 decomposes as {ip, A) with ip e r{T*B V2) and A e T{T*B Vq). The torsion of a generic 6* is 

T : B ^ V4 (8) (V: A V:) = V4 ® (V2 e Ve) 
(25) = (V2 © V4 © Ve) © (V2 © V4 © Ve © Vs © Vio), so 

T = (r| + r| + Tg^) + (r| + r| + t| + r| + rfo) . 

In this notation, each is a distinct irreducible component of T. The lower index n indicates 
the weight of the representation in which takes values, and the upper index r indicates the 
summand from which it was obtained. 

Thus, Cartan's first structure equation for a GL(2)-structure may be written in either vector 
form or polynomial form: 



(26) 



duj = -9 Alj + T{uj a uj) 
Explicitly, the connection term is 



(27) 



9 Auj^ 



f-8Lpo + X 8ip_2 

—2(p2 —4:(po + A 6(p-2 

—4:(f2 X 4:ip_2 

-6(p2 4v9o + A 

y 00 -8^2 



\ 




2v?-2 

i(/3o + Ay 



2 



A 



UJ 



so the matrix representation of 9 is (2(/?_2X — 2ipoH + 2ip)-£f -|- A/5), which takes values in ^1(2) C 
0KV4). 

Theorem 2.1. A generic GL{2) -structure tt : S — >■ M admits a two-dimensional family of connec- 
tions such that the essential torsion T decomposes irreducibly as 



(28) 



T = T2 + Te + Tg + Tio e V2 © Ve © Vg © Vio C V4 (8) (V4* A V4*). 



The two-dimensional family is parametrized by the possible equivariant inclusions 0/V2© Ve© Vg© 
Vio into V4 © (V4 A V4), and once such an inclusion is chosen, the connection is unique. 

Proof. Changes of pseudo-connection are of the form cp = + -P(t^) and A = A + Q{uj) where 
P e V2 © V| = V2 © V4 © Ve and Q e Vq © V| = V4. Preferred connections are obtained 
by analyzing the skewing map 5 that describes how changes of pseudo-connection (equivalently. 
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changes of horizontal section of the frame bundle) affect the torsion [1] : 



^ 0[(2)« ^ (V2 © Vo) ® VI ^ V4 ® (A^V:) — //°''(0K2)) ^ 0. 



(29) 




B 



For current purposes, Qi{2)^^'> and if°'^(0[(2)) are defined by the exactness of the sequence [T9] . 
They depend on the representation of the group GL{2), in this case the irreducible representation 
V4. To find the space of essential torsion, if°'^(0[(2)), one must compute SP and 6Q. 

FixP e V2©V4©V6, where V2 3 P{uj) = (P2, w)2 + (P4, w)3 + (i^6, c^)4- Let 5P G Va^A'^IV^) have 
components 6P = 6Pi + 5Pl + 5P| + ^Pa*' + 5Pf + 5P| + (JPg^ + 5Pfo, similar to the decomposition of 
T. Since these are irreducible components and the action of 5 must be 5'L(2)-equivariant, Schur's 
lemma implies 5 must preserve the weights of the representations. In particular there must exist 
constants 02, 04, ag, &25 ^4, and fog such that 5P2 = 02-^2, = ^2-^25 and so on. Thus, 

= {P{co),uj)^- 5P{uj,u) 

= ((P2, W)2 , W)i + ((A, C^)3 , W)i + ((^6, W)4 , ^)i 

- (a2P2, (w,w)3)o - (a4P4, (w,w)3)^ - (aePe, (w,w)3)2 

- (fe2P2, (c^, t<;)i)2 - {hPi, (w, 6^)1)3 - (fee^e, (w, w)^)^ . 

Carrying out this computation shows that 

(31) 02 = ^, 62 = ^, ^4 = ^, &4 = 0, ag = fee = 

Similarly, fix Q G V4, where Vq 3 Q{^^) = {Q'l,'^-')^- Let G V4 © A'^iV^) have components 
= 6Ql + 5Q| + 6Ql + + ^Ql + ^Qt + "^Qs + ^Qio^ but again the image must have the same 
weight as the domain. In particular there must exist constants C4, and such that 

= {Q{uj),uj)q - 6Q{u, u) 
= {{Q4, w)4 , w)g - (C4Q4, (w, uj)r^)^ - {d^Qi, {oj, w)i)3 . 

Carrying out this computation shows that 
(33) C4 = c/4 



40 160 

Fix a generic pseudo-connection {(p, A) with torsion T. Consider another pseudo-connection 
= + P{u}), A = A + Q{X) with torsion T. Then 

T(ci;, cj) = do; + ((^, u)^ + (^X, 

(3'^) = duj+{ip,u)^ + (P(w),cj)i + (g(cj),cj)o 

= (T + (5P + 5Q)(a;,a;). 
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Using Equation (31) and Equation (33), the absorption of torsion is dictated by the solvabihty of 
the following equations in terms of P and Q for fixed T and T: 

3 



(35) 



T 



n.2 
J- ft 



rpG 

^2 



To + — P2, 
2 10 ' 



-P4 



1 

40 



Q 



4, 



rpb 



rp6 

-^8 



rp6 



160 
20 



4, 



llQ — l^ 
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Generally, one may choose P2 to force exactly one linear combination of T| and T2 to vanish. 
Similarly, one may choose Pq to force exactly one linear combination of Tg and Tq to vanish. Unique 
Q4 and P4 eliminate T| and T|. All other components of T are fixed. Thus, Qi{2)^^^ = 51(2)'^^^ = 
and iJ0'2(s[(2)) = V2 © V4 © Vq © Vg © Vio, while H'^'\qI{2)) = V2 © Vg © Vg © Vio- □ 

Corollary 2.2 (Preferred GL(2) Connection). A GL (2) -structure B ^ M admits a unique con- 
nection A such that B has first structure equation 

dw = - (v3,w)^ - (A,a;)o 

(36) + (T2,(a;,a;)i)2 + (T6,(a;,a;)i)4 

+ (7^8, (t^, t^)i)5 + (T^io, (w, w)i)g . 



Proof. Of the possible connections, choose the one that absorbs T2 and T^. The remaining essential 
torsion is T = T| + T| + T| + Tf^. □ 



Henceforth, all references to 9, y). A, and T assume the connection in Corollary 2.2[ The specifi- 
cation of this connection over the others is arbitrary, but it does not affect any subsequent theorems 
in this article, since the V2 and Vg components of essential torsion turn out to be unimportant in 
the study of integrable PDEs. 



3. 2-lNTEGRABILITY 

Theorem 3.1. // a GL{2)- structure n : B ^ M is 2-integrable, then T = Tg (that is, T2 = Tg = 
^10 = Oj, and the bi-secant surfaces in M are parametrized by two functions of one variable. 
Conversely, if B is a smooth GL (2) -structure with T = Tg, then B is 2-integrable. 

Proof. To prove the theorem, one must find the conditions on B that allow an arbitrary bi-secant 
plane E G Gr2{TpM) to be extended to a bi-secant surface N C M. 

In a neighborhood M' of p, Fix u : M' ^ B{M'), a section of n : B ^ M; that is, fix u, a GL{2) 
coframe on M'. Let b = u{p). Since E is bi-secant, h{E) C V4 is spanned by (^^ix — hiyY and 
{g2X — h2yY with gih2 7^ (72^1- Through a GL{2) frame adaptation redefining m, one may assume 
that h{E) = span{a;^y^} C V4. Let E = u^{E) e Gr2(Tfei3). Then cj"^ A oj^\^ ^ 0. 
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Consider the linear Pfaffian exterior differential system X differentially generated by the 1-forms 
{cj~^, w^} with independence condition the Au^ ^ 0. It suffices to prove the existence of 
a surface N C B that is integral to X, because u{7i^{TN)) = a;„(TA^) = spanjx"^, y^} implies that 
the surface N = 7r{N) C M' is bi-secant. In fact, by adapting the moving coframe u appropriately, 
every bi-secant surface through E must arise this way. 

The generating 2-forms of X are 




(37) d|a;°| = | \ A {'^^^^, ] + \ \ uj~^ A oo^ 



modulo u ^,a;°,Ci;^, where 

r-2 = 48 T2,-2 + 8640 T6,_2 + 322560 T8,_2 - 4838400 Tio,-2 



(38) r'^ = -96 T2,o + 23040 Te^ - 4838400 T] 



10,0 



= 48 T2,2 + 8640 Tg^a - 322560 Tg^z - 4838400 Tio,2 

Because of the independence condition, integral elements exist only when the torsion can be 
absorbed. The torsion component r° can never be absorbed, so integral manifolds exist only when 
r° = 0. The condition of 2-integrability means that every 2-secant plane is tangent to a 2-secant 
surface, but the GL{2) action is transitive on 2-secant planes in TpM; therefore, 2-integrability 
implies that r° = for every element in the GL{2) orbit of T. Under the GL{2) action, the 
coordinates of the irreducible representations of T will change, so each irreducible representation 
that appears in r° must vanish identically. Hence, 2-integrability of M by integral manifolds 
implies 

(39) Tio = Te = T2 = 0. 

The remaining torsion components, r^^ and r^, are absorbed by setting tti = 2(y92 — 322560T8^_2 
and 712 = — 2v5-2 — 322560T8,2 i^""^, so 

(40) d I I = I I A ( ^ ) , mod tu 

This proves the torsion condition in the theorem. 

Conversely, to establish the existence and parametrization of bi-secant surfaces, one can apply 




Cartan's test for involutivity to the tableau in Equation (40) [51 For a generic flag of TpN 
obtained from generic linear combinations of w"^ and u'^, the tableau has Cartan characters si = 2 
and S2 = 0. Let 9 : N ^ B denote the embedding of the integral surface. The space of integral 
elements for the EDS (X, a;"'' Acu^) is 2-dimensional, as parametrized by the coefficients pi^4 and _4 
that appear in the pulled-back forms t'*(vri) = pi 4Z/*(a;'^) and i'*{7i2) = Ps ,-4^'*(i^~'^)- Therefore, if 
vr : i3 — 7- M is 2-integrable, then bi-secant surfaces in M depend on two functions of one variable. 
With the Cartan characters computed, Cartan's test for involutivity applies, so integral surfaces 
for the linear Pfaffian system exist in the real-analytic category. 

More can be said by employing modern theorems regarding hyperbolic exterior differential sys- 
tems [31]. The characteristic variety of the tableau consists of two real points, and because a 
generic tableau is involutive with si = 2, Yang's generalization of the Cartan-Kahler theorem to 
smooth hyperbolic systems implies that integral surfaces exist and are parametrized by two func- 
tions of one variable in the smooth category ^34t Theorem 1.19]. A special case of this observation 
is revisited in Section l64l □ 
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By restricting the torsion, Theorem |3.1| provides the necessary and sufficient first-order condi- 
tions for 2-integrabihty. These first-order conditions imply syzygies on the second-order invariants 
via the Bianchi identity. 

Corollary 3.2. Suppose the GL{2)- structure tt : B ^ M is 2-integrable with torsion T. Let S = 
V(T) denote the covariant derivative ofT, and let Q = T oT denote the symmetric product ofT. 
Then B has structure equations of the form 



dw = - (v3,w)i - (A,w)o + (T, 
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^4^^ ^'P = -\ + (^0. (^>^)3>o + (48Q4 + 42^4, {uj,uj),)^ 

+ (455*6, (w,w)i)5 + (SSS's, {uj,u})i)q + {-8Q4 - 12S'4, {uj,uj)^)^ 
dA = 960 {Se,{^,^)i)e 

for a scalar curvature function : B Vq- 

Proof. Write V = d + ^ for the covariant derivative on B defined by the connection 1-form A). 
Second-order consequences of 2-integrability arise from the Bianchi identity, V(6') Au = V(T{ijj A 
u)). Curvature appears in V(6'), which splits into R{ijJ Au) = dip + ^ {ip, ip)^ and r(a; A a;) = dA. 
The covariant derivative of the torsion two-form, V{T{ujAuj)), expands as V(T(ci;Aa;)) = V(T)(a;A 
(jj) + 2Q{T,T){u Au A u), so the Bianchi identity for a G'L(2)-structure is 

(42) R{uj A u) A u + r{uj A u) A u = V{T){uj A u) + 2Q{T,T){uj A u A u). 

Each of R, r, V(T) and Q is a function on B to the appropriate SL{2) module: 

R:B^si{2)®{VlAVl), 
r -B ^^®{yiAVl), 
^^^^ V(T) : B ^ i/°''(st(2)) ® VI, and 

Q:B-^ ^ym\H'^'\Q[{2)) n (V4 ® A'^VI). 



Since B is 2-integrable, T takes values only in Vg C if°'^(0[(2)). Using the Clebsch-Gordan 
decomposition, the irreducible components of these functions are (omitting the domain B for 
brevity) : 

R = {Rl + Rl + Rl) + {Rl + Rl + Rl) G V2 ® (V2 © Ve), 
r = r2 + r6 e Vo® (V2© Ve), 
^ ^ VT = ^4 + ^6 + ^8 + 5io + S12 e Vg © V4, and 

g = g4 + Qs G Sym2(V8) n (V4 ® (V2 © Ve)). 



Thus, Equation (42) and Schur's lemma together imply linear relations among the irreducible 
components listed in Equation (44). To find these relations, one can expand Equation (42) using 
the Clebsch-Gordan pairing; for example, one of the terms is 

(45) (VT, (w, w)^)^ = ((^4, w)o + {S&, + (^8, uj)^ + (5io, uj)^ + (^12, uj)^ , (w, . 

The resuh is ^10 = 0, i?^ = 33^8, Rl = 45^6, = 960^6, Rl = -I2S4 - 8Q4, Rj = 42^4 + 48g4, 
r2 = 0, and = 0. In particular, R^ is the only irreducible component of r and R that is not an 
algebraic function of T and V(T). □ 
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4. 3-Integrability 

Theorem 4.1. // an analytic GL{2)- structure rr : B ^ M is 3-integrable, then tri-secant 3-folds 
in M are locally parametrized by three functions of one variable. 



For 3-integrable (jL(2)-structures that arise from PDEs of hydrodynamic type as in Theorem 1.3 



(and are therefore also 2-integrable), this parametrization by three functions of one variable con- 
firms the computation presented in 



Proof. This theorem is proven by applying Cartan-Kahler theory to a differential ideal whose 
integral manifolds are tri-secant 3-folds N (Z M through an arbitrary tri-secant element E G 
Gr,{TpM). 

As in the proof of Theorem 3.1[ consider a local GL{2) coframe -u, and let h = u{p). Since E is 
tri-secant, b{E) is spanned by {gix — hiy)^, {g2X — h2yY^ and {g^x — h^yY, distinct, but u can be 
adapted so that b{E) is spanned by x^, y^, and (x + y^. Therefore 

(46) b{E) = {{A + B)x^ + B{4x^y + 6xV + W) + {B + C)y^ : A,B,C eR}c V4. 
Just as in Theorem 



3.1 



lift to E = u^{E) G G'r2(Tbi3), which is integral to the linear Pfaffian 
system X generated by /t^^ = — uo^ and = cj^ — uo*^ with the independence condition 
u:~^ A w*^ A 7^ 0. Again, the projection to M of any integral 3-fold of X will be a tri-secant 
3-fold that passes through E, and every tri-secant 3-fold arises this way (up to a GL{2) frame 
adaptation). 

The tableau and torsion of X are given by 



CO / a<b 



modulo K ^, K^, where tti = —2ip2, 112 = 2^5-2, and tts = 2y9_2 — 4(y9o + 4^92. The apparent torsion 
can be fully absorbed by redefining tTj = tTj — Pi^a^"' for the parameters 



(48) 



Pl,4 


- 

~ '-4,4 


P2-4 


- 

~ '-4,4 






P3,4 


~ '0,4 


P3-4 


= -Pl,-4 - 


P2,4 


= -P2,0 - 


Plfi 


= P3-4 - 



^2 _ 2 

"-4,4 ''"-4,0 



"0.4 = -P2fl + Tq! + T_l, + tL 



_ _ _ 2 _ 2 _ 2 

,4 ~ Pi -4 T-4A '^-4,0 '^"-4,4- 

The integral elements are still free up to arbitrary choice of three functions, pi-4, P2,o, and o- 
Since the Cartan characters are Si = 2, S2 = 1, and S3 = 0, but Si + 2s2 + Ss^ = 4 7^ 3, the tableau 
is not involutive; prolongation is required. 

Let X^^^ be the prolonged ideal, which is differentially generated by the forms k"^ and along 
with 

= TTi + pi _4 uj~^ - (pi _4 + r^4^4 + r^4_o + r^4_4) uj" - rr4_4 u^, 
(49) r]'^ = 712 + uj~^ + P2fl w° + (-P2,o + + Tlli + T^^) uj'^, and 
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After this prolongation, the tableau and torsion are given by 



(50) 





1 








\ 















1]^ = 




TXa 


-774 












7^5 


-775 


\rp j 


V 


-7r4 




) 



A 



modulo K^, T]^, rf, rf. This tableau has Cartan characters Si = 3, S2 = 0, and S3 = 0. Applying 
Cartan's test, Si + 2s2 + 853 = 3, which matches the dimension of the variety of three-dimensional 
integral elements of X'^^\ so the tableau is involutive. When the apparent torsion r*-^) vanishes 
or can be absorbed, then in the analytic category Cartan's test for involutivity implies that the 
integral 3-folds locally depend on three functions of one variable. 

In fact, one can say more by analyzing the characteristic variety of the involutive linear Pfaffian 
system presented in Equation (50). Given the highest non-zero Cartan character, si = 3, it is clear 
that the complex characteristic variety has dimension zero and degree three [3l Chapter V]. It is 
easy to verify that the characteristic variety consists of three real points, so each 3-fold integral to 
I^^^ is foliated by three families of 2-folds integral to A special case of this fact is revisited in 
Corollary |42 



□ 



To determine sufficient conditions for existence of tri-secant 3-folds, one must study the unab- 
sorbable portion of the remaining torsion, t^^\ Because the torsion of the prolonged system 

it will involve second-order invariants of the GL(2)-structure B that appear in the Bianchi 
identity for B, Equation (|42]). Since T is a priori valued in H°''^{gl(2)) = V2 © Ve © Vs © Vio, 
components of any of the following functions may occur in t^^^: 

R:B^V2® (a'V4), 



(51) 



r :i3->Vo©(A2V4), 
VT : i3 ^ (V2 © Ve © V; 



4, 



Vio) 

Q:B-^ Sym'(V2 © Ve © Vg © Vio) n (V4 ® A^Vi). 



The vanishing of the unabsorbable portion of r*^^) will place restrictions on the various irreducible 
representations appearing in Equation (51). The enormous complexity of Q and VT makes de- 
composition of the unabsorbable portion of t^^^ extremely difficult. Fortunately, one can make a 



simplifying assumption that is consistent with the motivating PDE theory in Theorem L3 hence- 
forth, all theorems discuss only those G'L(2)-structures that are both 2-integrable and 3-integrable. 

Corollary 4.2. If B ^ M is a 2, 3-integrable GL (2) -structure, then any tri-secant G M is 
triply foliated by bi-secant surfaces. Moreover, the net defined by TN (1 C is a coordinate net. 

The triple foliation by bi-secant surfaces appears to be new, though it is not surprising based 
on the description of highest-weight polynomial subspaces of V„ seen in [13]. This foliation shows 
that the "holonomic characteristic net" condition for tri-secant 3-folds seen in is superfluous 
since 2-integrability is implied for GL(2)-structures arising from Hessian hydrodynamic PDEs. 

Proof. Without loss of generality in a contractible neighborhood in M, one may apply a GL{2) 
change-of-frame so that TN = {Ax^ + B{x + y)^ -\- Cy^}. The characteristic net, TN fl C, is given 
by its tangent planes 

Ki = span{(a; + ^/)^^/^}, 
(52) K2 = spa,n{x^ , y"^} , and 

K3 = spanlx"^, (x + y^}- 
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Each of Ki, K2, and K3 is clearly bi-secant. It suffices to prove that they are everywhere tangent 
to the level sets of a coordinate system on A^, as these level sets will be bi-secant surfaces. 

Recall that N is define d by the projection of an integral manifold of the prolonged system X^^^ 



on i5 X M'^ from Theorem . In particular, there is a submanifold N G B with embedding map 
V : N ^ B depending on three parameters Pi, P25 -P3 such that each TpN is the projection of 

TfeiV = ker ({/^-^fi:^r/^r/^^7^ 
7r4-Pi(a;"^-w°), 
7r5-P2(c^°-u;'), 
7r6-Pia;-^ + P3a;°}). 



(53) 



Let = ^ — = w^, and = cj^ — . Therefore, at each basepoint, Ki C TpN is the 

projection of the plane ker(i/*(T*)) C T^N. One can easily compute that d(z>*(T*)) = modulo 
i/*(T*) for each i; therefore, the Frobenius theorem provides local coordinates (s^, s^, s^) on such 
that Ki = ker(ds^). □ 

When B arises from a Hessian hydrodynamic PDE, these coordinates are called the Riemann 
invariants of the hydrodynamic reduction. 

Theorem 4.3 (2,3-integrable (jL(2)-structure equations). A GL (2) -structure B is 2,3-integrable 
if and only if the torsion T of B only takes values in Vs and the curvature is a function of T . In 
particular, every 2,3-integrable GL{2)- structure B has the following structure equations 

du = - {(f, uj)-^ - (A, lu)q + (T, (w, u)^)^ 
dA = 

= ^)i - 2080 ((T, T)s , {u, u),)^ + 64 ((T, T)^ , {u, u),)^ 

^^^^ - y ((^> T)6 , (w, w)i)4 + y {{T, T)^ , (w, a;)^)g 

dT = J{T) A 

for a 9 X 9 matrix-valued function J whose entries are linear and quadratic polynomials in the 
coefficients ofT, as provided in Appendix\A\ and 



Proof. 2-integr ability implies that T = Tg, so P, r, V(T), and Q decompose as in Equation (44) 



with the relations implied in Theorem 3.1 3-integrability implies the vanishing of the GL{2) 



orbit of the unabsorbable portion of r(i) in Equation fco|. The vanishing of r^^^ and the Bianchi 
identity (d^ = 0) together imply the additional relations S12 = 0, 75*8 = 24(^8 = 24 (T, T)^, 5*6 = 0, 
21^4 = = 8 (T, T)g, and P^ = -2080 (T, T)g. 

Therefore, the curvatures, P and r, are quadratic functions of T. The derivative of torsion, dT, 
is also a quadratic function of T, as expressed in the matrix J{T). Hence, T is the only invariant 
of any order for 2,3-integrable G'L(2)-structures. □ 

Definition 4.1. The notation {B, M,p)2,3 indicates a smooth 2,3-integrable GL{2)- structure it : 
B ^ M such that p & M and such that M is connected. 



The condition that M is connected is crucial in what follows, and pointedness is technically 
useful. 
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Definition 4.2 (Representatives). {B, M,p)2,3 represents f G Vs if v & T{Bp). More generally, 
a 2,3-integrable GL (2) -structure B ^ M represents v E Vg if T{B) 3 v. Likewise, a Hessian 
hydrodynamic PDF represents v E Vs if the induced 2,3-integrable GL (2) -structure over M = 
F~^{0) represents v. 

Lemma 4.4. The singular distribution on Vs defined by the columns of the matrix J is integrable, 
providing a stratification of Vg into leaves that are submanifolds . That is, for any v G Vg, there 
exists a unique submanifold Oj{v) such that T„Cj(t>) = range J(f). 



Lemma 4.4 is elementary in its modern interpretation using the theory of Lie algebroids and 



smooth groupoids. Since the equations in Theorem |4.3| are closed under exterior derivative, the 
matrix J is closed under the corresponding bracket, so it defines the anchor map of a Lie algebroid 
over Vs. This Lie algebroid is neither regular nor transitive, but there exists an integrating groupoid 
over the base. Vs. This groupoid is smooth, but it may not be a Lie groupoid. However, the 
groupoid is transitive when restricted to each of its orbits, Oj{v), and each orbit is a submanifold 
of the base [121 1231 [2E]- These orbits may be regarded as the leaves of a singular foliation of Vs 
such that TyOj{v) is spanned by the columns of J{v). 

Corollary 4.5. For any v G Vs, there exists a real-analytic connected 2,3-integrable GL{2)- 
structure B ^ M such that T{B) 3 v. That is, every v E Vg is represented by some real-analytic 
{B, M,p)2^3. Moreover, J has constant rank on any such B. 



Proof. Since the structure equations in Theorem 4.3 are closed under exterior derivative, this is a 



direct application of the existence part of Cartan's structure theorem, which is an generalization 
of Lie's third fundamental theorem to the intransitive case [S]. See Appendix A of [6] for a clear 
summary of the special case of Cartan's structure theorem needed here. As used here, Cartan's 
structure theorem holds in the smooth category; however, because T o T and J{T) are real- 
analytic (in fact, linear and quadratic), one can use the Cartan-Kahler machinery to produce a 
local real- analytic solution B of the structure equations. On any solution B, the image of dT\b is 
the image of the columns of the matrix J{T{b)), and these columns span the tangent space of the 
submanifold Oj{T{b)). Thus, on any connected B containing b, the rank of J equals the dimension 
of Oj{T{b)). □ 

Corollary 4.6. (S, M,p)2,3 and (i3,M, p)2,3 admit a local GL{2)- equivalence f : M ^ M with 
fijp) = p if and only ifT{Bp) r\T{Bp) ^ 0. That is, the value T{b) uniquely defines (i5, M, p)2,3 in 
a neighborhood of p = 7r(6), and this local GL{2)- structure is real- analytic. 

Proof. This is a direct application of the local uniqueness part of Cartan's structure theorem. 
Again, the theorem applies here in the smooth category, but the real-analyticity of the structure 
equations implies that any 2,3-integrable G'L(2)-structure is locally equivalent to the real-analytic 



GL(2)-structure produced in Corollary 4.5 □ 



Corollaries 4.5 and 4.6 establish the primacy of T in the study 2,3-integrable GL(2)-structures 



and the related PDEs. However, Corollary 4.6 provides only local GL(2)-equivalence, so it is useful 
to have a weaker "chain-wise" notion of equivalence that applies to global (but still connected) 
GL(2)-structures. 

Definition 4.3 (Leaf-equivalence). {Bq, Mo,po)2,3 O'lT'd {Bk, Mk,Pk)2,3 o,re said to be leaf- equivalent 

if there exist finite sequences {(i3j, Mj,pi)2,3} and [vi] with l<i < k—1 such that (i3j, Mj,pj)2,3 
and {Bi+i, Mj+i,pj+i)2,3 both represent Vi for < i < k—1. 



The term "leaf-equivalence" arises from the leaves of the singular foliation of Vs from Lemma 4.4 



These leaves separate all possible (i3, M,p)2,3's into equivalence classes by the value of T. By slight 
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abuse of notation, if (i3,M, p)2,3 represents v, then write Oj{B) to denote Oj{y). Leaf-equivalence 
can now be rewritten more succinctly. 

Theorem 4.7 (Leaf-equivalence). (S, M, p)2,3 and {B, M,p)2,3 are leaf- equivalent if and only if 
Oj{B) = Oj{B). Moreover, T : B Oj{B) is a submersion. 

Although T : B ^ O j{B) is an open map, it need not be surjective. 



5. The Classification 

To classify connected 2,3-integrable (j>'L(2)-structures is to explicitly identify the leaves of the 
foliation of Vs. To identify the leaves of the foliation is to identify connected GL(2)-invariant sub- 



manifolds of Vs where the rank of J{v) is constant. This is all ultimately achieved in Theorem 5.3 
thanks to the observations in Lemmas 15.11 and 15.21 

Lemma 5.1. In vector space of l^th degree polynomials in the variables f_8, . . . ,f8, the determi- 
nant of J{y) is a non-zero scalar multiple of the discriminant of the polynomial v G Vs. Therefore, 
J{v) is non-singular if and only if v has eight distinct roots. Moreover, if v is a nontrivial polyno- 
mial with k distinct roots, then the rank of J{v) is k-\-l. 

Proof. The determinant statement is verified by computing the g.c.d. of the two polynomials, 
det J{v) and disc(f). The rank statement is directly verified by writing v = {gix — hiy){g2X — 
h2y) ■ ■ ■ {gsx — hsy), imposing multiplicity on the /ij's and gi^s and computing the rank of J{v) 
directly. □ 

As it happens, a simple exercise in linear fractional transformations shows that the multiplicity 
of the roots is preserved by the GL{2) action on Vs [H ITT] . 

Lemma 5.2. For v{x,y) G Vs, the multiplicity and complex type of the roots are preserved under 
the irreducible action of GL{2,'R) . 

Bearing in mind this lemma, notational shorthand for the root type of a polynomial v{x, y) G Vs 
is useful: Suppose v{x, y) factors as {gix—hiyY^ ■ ■ ■ {gmX—hmyY"' such that ri > r2 > ■ ■ ■ > and 
ri-\-r2+- ■ ■+rm = 8. These exponents define a partition of 8 that is written as {ri, r2, . . . , Vm}- If the 
roots gk/hk and gk+i/h^j^i are complex conjugates, then denote this by enclosing their exponents 
in square-braces: {ri, . . . , [r^, r^+i], . . . , r^}. Denote the root type containing v{x, y) by [v{x, y)]. 
A root type is equivalent to the subset of Vs comprised of all polynomials that factor according 
to the given partition and complex-conjugate pairing. For example, [x'^{x + iy^lx — iy^] = 
[{x — 5yY{x + 2iyY{x — 2iyY] = {4, [2,2]} C Vs. Let {0} denote the trivial root type, the zero 
polynomial. There are 54 non-trivial root types that partition Vs into strata from dimension two 
to dimension nine, as represented in Figure [Tj 

In Figure [T| arrows mean "closure contains." Shaded nodes represent root types that contain 
exactly one GL{2) orbit. Oval nodes represent open root types. The square node represents 
the nearly-closed root type, {8}, which (when is included) is the rational normal cone in Vs. 
Hexagonal nodes represent root classes that are neither closed nor open. Note that strictly real 
root types actually have two connected components (for example, {8} is comprised of the two ends 
of the rational normal cone), but this is another artifact of the projective nature of the group 
action, as the two ends are in the same GL{2) orbit by the — Jg action. 

Theorem 5.3 (Leaf-Equivalence Classes). The leaf- equivalence classes of connected 2,3-integrable 
GL (2) -structures are exactly the root types in Vs. That is, for all v G Vs, Oj{v) = [v]. 
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Figure 1. The stratification of Vg into the 55 root types, sorted by dimension. 
Shaded nodes indicate root types that contain exactly one GL{2) orbit. The square 
node is closed (if is included); the oval nodes are open, and the hexagonal nodes 
are neither closed nor open. 



Proof. For each f , both Oj{v) and [v\ are smooth submanifolds of Vs, and by Lemma 5.1 they 
have the same dimension. To prove that Oj{v) = [v], it suffices to prove that T^Oj{v) = T^[v] 
for all V eVs- 

Fix a root type [v] and an arbitrary f G [v]. Note that 



T.OjiB) = DT,{T,B) = range J(y). 



Thus it suffices to find, for each column Ji{v), a tangent vector D{v) G T^[v] such that D{v) = 
Jiiv). As it happens, the equations defining the various D{v) are easily solvable at arbitrary points 
in all 54 non-trivial root types and for all columns of J. Because both the leaf-equivalence classes 
and the root types partition Vg by smooth submanifolds, and because T^Ojiv) = Ty[v] for all 
w G Vs, the partitions must be identical. 



20 



A. D. SMITH 



To illustrate the computations, consider an arbitrary point v 
[v] = {8}. An arbitrary element of T^[v] looks like 



{gx + hyY in the root type 



(55) 



D{v) = 8{Gx + Hy){gx + hy)"' 



( ^Gg' \ 

Hg'^ + 7Ghg^ 
2Hhg^ + GGh'^g^ 
SHh^g^ + 5Gg^h^ 
AGg^h^ + 4Hg^h^ 
5Hg^h^ + 3Gh^g^ 
GHh^g^ + 2Gh^g 

IHh^g + GK^ 

\ mh' J 



Therefore, one must solve Div) = JAv) for G and H in each of the columns i 



, 9. In this 



case, 



(56) 



J{v) 

















g' 


-16V 


16^^ 



















hg'' 


-14/iV 




2g^ 

















h^g"^ 


-12h^g^ 
-10/iV 


8/^2(^6 


Ahg'^ 

















h^g^ 


Ah^g^ 


Gh^g'' 

















h^g^ 


-8/iV 





8h^g^ 

















h'g' 


-6h^g^ 


-Ah^g^ 


10/tV 

















h^g^ 


-AK^g 


-8/iV 


12h^gi^ 

















h'g 


-2h^ 


-Uh'^g 


lAh^g"^ 


Vo 



















-16h^ 


Wh'^g / 



For each column Jj(f), it is easy to see the solution values of G and H. 

For the other root types, the computations are similar but somewhat more complicated. All 
that matters is the fact that they can be solved for arbitrary v. □ 



5.1. Symmetry Reduction. Geometrically, it is interesting to reduce all symmetry from a struc- 
ture. Fix {B, M,p)2,3 and v e T{B). Let B'' = {b e B : T{h) = v}. This is a sub-bundle of B with 
fiber group Stab(f) C GL{2). If dim[f] < 3, then the stabilizer group is smooth; however, the 
stabilizer groups must be discrete for all the larger root types. The stabilizer group, either smooth 
or discrete, must always appear in the well-known list of GL(2, C)-stabilizers of polynomials [2]. 
In the discrete case, the fiber must therefore be either a cychc group or a dihedral group [27] . 

The structure equations for B"" show no dependence on T (as it has been fixed), so a neighborhood 
in B^ is a local Lie group of dimension 9 — dim[f], and M is locally the homogeneous space 
;B''/Stab(f ). The relations defined by kerdT(f) = ker J(f) determine the pull-backs of u, A, and 
if to B"", so one may explicitly reduce Equation to obtain structure equations for the local Lie 
group B". 

If f = 0, then Stab(f) = GL{2), so the local flat 2,3-integrable GL(2)-structure is the local Lie 
group of dimension 9 obtained by setting T = in Equation (54). 

The root type {8} is a single GL{2) orbit, so any representative will generate all representatives. 
Suppose V = x^. Then dT_8{x^) = A + IGipo and dT_e(x^) = 2(f2 while dT^ = for /c > —6. Since 
dT is a vertical 1-form on B, the value of T only varies in the fiber of B. The reduced structure 



GL(2) PDES 



21 



equations are 

dw-^ = 24 (/?o A uj-^ - 8 (/?_2 A + 2 • 322560 A w^ 
dw"2 = 20 (/^o A uj-"^ - 6 A + 322560 cj^ A w^ 

(57) V- ^ ' 

dw^ = 12 A - 2 V9_2 A cj^ 

4 o ,„ A , ,4 



dw = 8 A u; , 
d^^o = 0, dv9_2 = 4 A 

These equations can be easily integrated, so B^^ has coordinates a, and 6 such 

that 





= a ^ da 




= db, 


UJ^ 


= a' de^ 




= a'' {de 




= a'^ {de 




= a'' (dr 




= a''{dr 



(58) ^ ^ 

uj" = a" [^ai^' — 4 

66 d^ + 126^ d^ - 8b^ d^ - 322560^^ d^^) 

36 dr^ + 246^ d^ - 32b^ d^ + 166^ de^+ 

8(322560)^^6 d^^ - 2(322560)^^ d^). 

The field of rational normal cones can now be written in these local coordinates, since Cp = 
{{up{v) , Up{v)) 2 : V G TpM} for a local section u of B. 

Each of the 3-dimensional root types is a single GL(2)-orbit. Again, T{B) = [v] for any v in 
the root type, so an arbitrary representative v may be chosen for any {B, M,p)2,3. Each v has 
a 1-dimensional stabilizer, which is a Lie subgroup of GL{2). The corresponding Lie algebras 
are easy to compute |2l[27j. In all cases, dT has both vertical and semi-basic components, so the 
embedding of the stabilizer fiber group varies over M. The reduced structure B"" is a 6-dimensional 
Lie group. 

If dim[f ] = 4, then B"" is a finite cover of M, but there are only eight such structures, since each 
root type with three roots is a single GL(2)-orbit. This contrasts with the case dim[f] > 5. If 
dim[f ] > 5, then B"" is a finite cover of M, but there is no reason to believe that T{B) C [v] implies 
T{B) = [v] when dim[w] > 5. Consider [v] = {2,2,2,2}. The quotient space {2, 2, 2, 2}/GL(2) has 
orbifold singularities; for example, the point x^y'^{x + y)^(a; — yY has an eight-element stabilizer 
group, but all nearby points have trivial stabilizer groups. The existence of these orbifold singu- 
larities implies that there cannot be a surjective smooth map M — )■ [v\/GL[2). If T{B) is a proper 



subset of [v] for all B representing f , then by Theorem 5.3 a finite sequence of 2,3-integrable GL{2)- 
structures connects any two points in the leaf, but infinitely many locally distinct G'L(2)-structures 
are required to cover the entire leaf. This behavior is closely related to the topology of CP" with 
k marked points, up to GL{2, C) action, which is a difficult and well-known problem in complex 
algebraic geometry. The real case encountered here is both harder and less well-known than the 
classical complex case. The most relevant result is [L0\, which shows that the orbifold singularities 
present in the larger leaves of Vs are so bad that they preclude the existence of Riemannian metrics 
on the leaves. 

In the open case, dim[f ] = 9, so J(T) has maximum rank, and T itself provides local coordinates 
on B. 

This discussion is summarized by the shape of the nodes in Figure [2] 
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6. GL(2) PDEs 

This section contains various conclusions regarding Hessian hydrodynamic PDEs that can be 
inferred from Theorem L3 and Theorem 5^ Recall that a PDE F = is said to have k symmetries 
if the Lie algebra of point symmetries of F~^{0) has dimension k. 

Lemma 6.1. A Hessian hydrodynamic PDE F{uij) = representing f G Vg has k symmetries if 
and only if dim[f ] = 9 — k. 

Proof. The contact transformations by CSp{3) induce the local automorphisms on B, but an 
automorphism ip : B ^ B near 6 is a symmetry if and only if the structure equations are preserved 
by ip*. This can happen if and only if T{b) is preserved near ip{b). In particular, ilj*{dT) = dT if 
and only if tp* is the identity on the range of dT, which has dimension dim[f]. □ 

Corollary 6.2. There is no Hessian hydrodynamic PDE with exactly eight symmetries. 

Proof. There is no root type of dimension one. □ 

The classification provides even more bountiful information about the Hessian hydrodynamic 



PDEs, because Theorem 4.3 allows a converse of Theorem 1.3 



(59) = _L , tt = a*, 7 = 7*, d/i + /i A /i = 0. 



Theorem 6.3. Every {B, M,p)2,3 is realized, locally nearp, by a Hessian hydrodynamic PDE. 

Proof. The structure {B, M,p)2^3 is realized, locally near p, by a Hessian hydrodynamic PDE 
as in Lemma 1.2 if and only if there exists a neighborhood M' (Z M oi p and an embedding 
i : M' — )■ A° such that the distribution of rational normal cones z=k(C(M')) over i{M') is the same 
as the intersection of i*(TM') with the distribution of Veronese cones over A°. 

Since the fibers B are exactly the symmetries of C and the fibers of CSp{3)° are exactly the 
symmetries of the distribution of Veronese cones, it suffices to establish a bundle immersion h : 
B{M') CSp{3y covering an embedding i : M' ^ A°. 

Let n denote the Maurer-Cartan form of CSp{3), which is of the form 

7 

(One may assume that the conformal scaling has been incorporated into /3, as it is below.) By 
Equations (16) and (17), the components /3 and 7 are vertical for H, and a is semi-basic. 

Recall the Fundamental Lemma of Lie Groups [22l Theorem 1.6.10]: If there exists rj : TB — >■ 
csp(3) such that d?7 + 77 A 77 = 0, then for any b & B, there exists a neighborhood B' of b and a 
map h : B' ^ CSp{3) such that h*{fi) = rj. Moreover, if h*{a) is semi-basic, then the fibers of B' 
immerse into the fibers of CSp{3). 

Therefore, it suffices to construct a C5p(3)-valued Maurer-Cartan form rj on B such that the 
entries of the lower-left symmetric submatrix of rj, namely h*{aij), are semi-basic on B. For 
brevity, the pull-backs h*{-) are dropped from the notation henceforth. 

The condition that a is semi-basic is that aij = AijaOj"". Note that a : V4 — )■ Sym^ 
Sym^(V2), and recall that for u, f G V2 the symmetric tensor is given by 

Therefore, to respect the weights in the GL{2) representation, a must have the following form for 
constants A_4, A_2, Aq, A'q, A2, and A4: 

(61) a = {A^2 OJ-^ A'o w° ^2 

Aq 00° Ao 00'^ A4 00^ 
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Writing (3 = (3y,{ip) + f3x{X) + (3t{^), it is apparent that (3^ must be the representation sl{2) — >■ 
M3x3(M) such that the natural action of M3x3(M) on Sym^(V2) is induced by the natural action of 
s[(2) on w G V4. One can now easily verify the following formulas: 



a 

(62) 

(3 = I -4v9_2 4(^2 I - ^Xh + PAoj) 

Here, the — IA/3 component of /3 is simply the scaling action of GL{2) as represented by the scaling 
action in CSp{3). If PGL{2) and S'p(3) were used instead, it would not appear. 

All that remains is to find (3t and 7 such that dr] + r] /\ r] = Q . This is arithmetic, and solutions 
exist. The simplest 77 (the one where the undetermined coefficients in Pt are set to 0) is provided 
in EEI. □ 





00-' 




















2V92 


-4v9- 












-2y^_2 



Theorem |6 . 3 1 allows explicit construction of several of the most-symmetric Hessian hydrodynamic 
PDEs. Each of the root types with three or fewer roots is itself closed under GL{2)\ hence, the 
torsion of any G'L(2)-structure covers the entire root type. This allows a first step at constructing 



all Hessian hydrodynamic PDEs. However, Theorem 6.3 does not say that there are only 55 
Hessian hydrodynamic PDEs up to (75*^(3) actions. There is no reason to believe that the torsion 
map T : i3 — > Vs is surjective for the larger root types. When T fails to be surjective, there cannot 
be a single PDE that represents all possible torsions in its root type. Nonetheless, the leaves of 
dimension at most 4 can be used to fully describe the most symmetric Hessian hydrodynamic 
PDEs, as in the following theorems. 

Theorem 6.4. The root type {0} is represented by the wave equation, 

(63) U22 = uis 

and this representation is unique up to CSp{3). Therefore, the wave equation is the unique Hessian 
hydrodynamic PDE with nine symmetries. 

Theorem 6.5. The root type {8} is represented by the first flow of the dKP hierarchy, 

(64) U22 = Mi3 - ^ (^33)^ 

and this representation is unique up to CSp{3). Therefore, the first flow of the dKP hierarchy is 
the unique hydrodynamic PDE with seven symmetries. 

Theorem 6.6. The root type {7, 1} is represented by the PDE 

. , 11 

(65) U22 = ui3 - — M33 + -M33M23 

and this representation is unique up to CSp{3). 

Theorem 6.7. The root type {6, 2} is represented by the PDE 

(66) U22 = ui3 + — 

5M33 - 14 

and this representation is unique up to CSp{3). 
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Theorem 6.8. The root type {6, 1, 1} is represented by the PDF 
, 7m23 (m23 - M33) , 49(-(m33)^ + 14m33 - 28) 

(67) ti22 = Ui3 H 7^ 7-;^ h 



(5^33 - 14) 12(5^33 - 14) 

and this representation is unique up to CSp{3). 
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Theorems 6.4 through 6.8 are proven via the same technique. For illustration, consider the 



simplest non-trivial case, Theorem 6.5 



Proof of Theorem\6.5. Consider a 2,3-integrable GL(2)-structure vr : i3 — )■ M with T{h) 



X 



1 



322560' 



SO [r(i3)] = {8}. The constant is chosen for the aesthetic appeal of the resulting PDE. The goal is 
to describe Theorem 6.3 s embedded hypersurface i{M') = {Il{h{q)) G A° : g near b} as the locus 
of a single equation F{U) = 0. 

As CSp{3)° is a matrix group, the open set {g = h{q) : q near b} can be described as {expj{rih{v)) : 
V G TbB}. Moreover, only v G ker(y9. A) need be considered, since the fibers of vr : S — M immerse 
into the fibers of 11 : CSp{3)° — ?■ A°. Fix an arbitrary v G ker(y9. A) C TbB and write v in compo- 
nents using the tautological 1-form, v^u = (f_4, f _2, fo, i'2, ^^4) G Of course, v^u does not 
actually provide local coordinates on B or M; however, the matrix rif){v) still represents a generic 
point in h^(Ti,{B)), as seen here: 



(68) 



Therefore, 



Vb{v) 



( ° 































































V4 


V-2 


Vo 


V2 











\ Vq 


V2 


V4 








0/ 



/ 



(69) 



1 


-V4, 



v-i - Ivl 

V^2 - \V^V2 
1 2 
^^0 - 2^4 



V-2 















o\ 


1 

















1 











+ ^f4f2 


Vo + \v\ 


1 









V2 





1 





V2 













So, using Equation (17), a generic point U G i(M') looks like 

V4 



6^4 



+ (^^0 + \VI)V^ V_2 + \V\V2 



(70) 



U = U{expj{r]p{v))) 



1, 



f_2 + ^V4V2 

Vo + \vl 



Vo 
V2 




(71) 



^722 = U- 
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33; 



There is a single relation between the entries of such U: 

1 

~ 2 

When U is interpreted as the Hessian of u : M'^ — t- M, Equation ( 71 ) is the first flow of the dKP 
hierarchy, a well-known example of a Hessian hydrodynamic equation. □ 

The only change for the other root types is that 77 is more complicated; hence, its exponential 
is (immensely) more difficult to compute, and the relation F{U) = is more difficult to recognize. 
Note also that f/ = is always in the locus of the equation obtained by this procedure. Therefore, 
PDEs such as the Boyer-Finley equation, u^x + Uyy = e''" (which has six symmetries and must 
represent one of the 3-dimensional root types [E]), will not directly appear as representatives 
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via this procedure. None-the-less, every CSp{3) equivalence class of Hessian hydro dynamic PDEs 
must arise this way. 

6.1. Hyperbolic Planar PDEs. This section presents some preliminary but intriguing observa- 



tions regarding the hyperbolic linear Pfaffian system I describing bi-secant surfaces in Theorem 3J^ 
and its relation to hyperbolic second-order planar PDEs, 

(72) /(^\ ^^ z, zi, Z2, Zn, zi2, Z22) = 0. 



Equation (72) defines a 7-dimensional manifold Ej = /~^(0) C J^(M^, M) whose structure equations 
are obtained by pulling back the contact system [181 [30]. Such admit a point- wise classification 
into Monge-Ampere, Goursat, or generic-type equations. 

Theorem 6.9. Consider {B, M,p)2,3 with T{b) = v for some b G Bp. Over a neighborhood M' of 
p, there is a bundle W — )■ M' with 7-dimensional fiber and a submersion f : B{M') — )■ W such 



that the ideal X from Theorem 3.1 describing the existence of bi-secant surfaces through p is the 
pull-back of a hyperbolic linear Pfaffian ideal X on W . This W admits a coframe {P^, . . . , /?'') such 
that 

d/3^ = AP^ + P^ A P\ mod P^ 
(73) d/32 = Ui P^ AP'^ + P^ A /3^ mod P^ 

d/3^ = U2 P"^ AP^ + P^ A P\ mod P\ p\ 

for Ui = — 645120T_8^^^ and U2 = — 645120T8^^^ for some non-zero functions /ii, fi2, and fi^ 
on W. 

Moreover, there is some second-order hyperbolic planar PDF f and a diffeomorphism : — >■ 
W such that these structure equations on W pull back via ^p* to the contact-induced structure 
equations on J]f. 



Proof. Fix {B, M, p)2,3- As in the general case of Theorem 3.1 , the linear Pfaffian system describing 
the existence of bi-secant surfaces through p is differentially generated by o;"^, a;°, and a;^, and it 
has tableau given by 

,-2^ 



(74) 





for TTi = 2(p2 — 322560T_2Ci;^ and 112 = — 2(/?_2 ~ 322560T2a;~^. The Lie algebra A{X) of Cauchy 
characteristics of X is spanned by the duals of A and ip^. Therefore, the retracting space C(X) = 
A{X)-^ is a rank-seven Frobenius system on B, so B admits a foliation by 2-dimensional Cauchy 
characteristic surfaces [3l Section 11.2] [221 Section 6.1]. Let W denote the 7-dimensional (local) 
leaf space for this foliation, so there is a submersion fc : B ^ W. Since the u is semi-basic for the 
submersion vr, W also admits a submersion onto a neighborhood of p G M. It remains to find the 
structure equations for a coframing on W. 

Write = 0;°, = o;"^, = w^, = tti, = u'"^, = 712, a'' = OJ^, = fo, and 
= A as a coframe for B. To simplify the notation, fix the index convention 1 < i, j, k < 7 and 
8 < r, s < 9, so the forms a* are semi-basic for the bundle B ^ W, and the forms are vertical 
for the bundle B -^W. Write da* = -^Ci^ A a'' - C'j.^.a^ A a*", where Cj;, = -Clj and Cj„ are 



functions of T as determined by Theorem 4.3 Note that Cj^ = if i 7^ j, and 

(75) Cl, = 0, C'2s = -4, = 4, Cts = 4, Cl, = -8, Cl, = -4, C',, = 



(76) 
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Corollary 2.3 of |3] implies that there exist functions /ij on B such that —a" 



^J.i 



is basic. Define 



a new co-framing (a*) for B by setting a* = —a* (no sum) and = . Thus, T*i3 is (locally) 
split into basic 1-forms and vertical 1-forms with respect to the bundle vf : — ?■ W . Of course, 
a* = 7f*(/5*) for some independent 1-forms and (/3^,. . . ,/?'') is the desired co-framing for W . 
Since dt is basic, the structure equations of (a*) and (/?*) are identical. 

There is a lot of freedom in the designation of /ij. Write d/ij = + fii^rC^"^ . The condition that 
(5* is basic is equivalent to fii^r = f^iClr for all 1 < i < 7 and 8 < r < 9, but /ijj- is otherwise free, 
so any non-zero function /j on W may be lifted to B in local coordinates (x^, . . . ,x^) by setting 
/^i = /i Gxp(C'iga;* -|- ClgX^). 

At this point, the coframe (/?*) satisfies 



(77) 



d/3^ 


= 0, 


d/32 








d/33 









mod {I3\ 13^/3'} 



so one may effectively eliminate /X5 and /i7 by redefining as 
structure equations for W satisfy 

d/3^ = 0, 

(78) d/32 = A /3l 

d/3=^ = A 



and as 



Thus, the 



mod {/3^/3^/33} 



One may now re-label the coframe (/?*) following the procedure given in Appendix A of [3D] to 
obtain Equation (73). 

Although Ui and U2 as written in the theorem are not explicitly functions on W, this is easily 
remedied. As noted above, Hi may be taken as the lift of fi on W. Also, the action of and 
is diagonal on T: 



(79) 



/T-s\ 


















12T_6 


T_4 




T_4 




8T_4 


T_2 








4T_2 


To 




To 


a® + 





T2 




T2 




-4T2 


T4 




T4 




-8T4 






Te 




-12T6 


\Ts/ 




\TsJ 







mod a^, 



, a 



So, while there is no natural map from W to Oj{B), the local Lie group G generated by the fiber 
actions of H and Jg (corresponding to and a^) induces a map W — )■ Oj{B)/G that one could 
also call T. 

Theorem 11.1.1 of |29] implies that any 7-dimensional manifold with structure equations of this 
form must be diffeomorphic to for some hyperbolic planar PDE /. □ 

Corollary 6.10. For any Hessian hydrodynamic PDE, the level sets of the Riemann invariants 
are the solutions of planar hyperbolic PDEs. 



Proof. The manifold W arises from the ideal describing bi-secant surfaces in M, and Corollary 4.2 



shows that these are the level sets of the Riemann invariants of the Hessian hydrodynamic PDE 
defining B. □ 
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This corollary is not at all surprising, since the entire point of hydro dynamic reduction is to 
reduce a hyperbolic PDE in three variables to a family of hyperbolic planar PDEs defined by the 
Riemann invariants. 



In the case of 2,3-integrability, Corollary 6.10 also provides a more explicit justification for the 



claim in Theorem 3.1 that the PDE defining bi-secant surfaces can be solved in the smooth category 



with smooth initial data. One naturally asks "to which T,f is W equivalent?" 

Corollary 6.11. Let T,f = f^^{0) C J^(M^,]R) for a hyperbolic planar PDE f. Suppose there is a 
local diffeomorphism ip : T,f ^ W with Tx{h) G ipijlf). 

(1) If X = and y = are both roots of T{b), then is of the Monge-Ampere type at 

r\m), 

(2) // exactly one of x = or y = is a root of T(b), then is of the Goursat type at 

(3) // neither x = nor y = is a root of T{b), then is of the generic type at ip ^(7r(6)). 
Proof. These are the three pointwise types of planar hyperbolic PDEs, and for structure equations 



of the form in Equation (73), they are determined by whether Ui and U2 vanish [TRi ISO] . Since 
Ui ~ T_8, f/i = if and only if y = is a root of T{b) G Vs. Similarly, f/2 = if and only if x = 
is a root of T{b) G Vs. □ 



Corollary 6.11 puts interesting restrictions on which GL(2)-structures can yield which planar 
PDEs. In particular, x = and y = Q can both be roots of f = T{b) if and only if [v\ is a root type 
having two distinct real roots, and the strictly complex root types cannot have x = or ?/ = as 
roots. Notably, the type of W can change. For example, suppose T{b) = x'^y, so '?/'~^(7r(6)) is of the 
Monge-Ampere type but nearby ip'^^Ti^b')) is of the generic type, since T{b') = {x + siyY {e2X + y) . 
This type-changing does not occur for the flat structure, and it is easy to compute a change-of- 
frame from the structure equations of the fiat W to the structure equations of for the planar 
wave equation 2:12 = 0. 

Corollary 6.12. Let {B,M,p)2,3 be fiat. Then W is isomorphic to {zi2 = 0} C J2(]R^]R). 

7. Concluding Remarks 

The main results of this article are summarized in Figure [2j In short, Hessian hydrodynamic 
PDEs in three independent variables are equivalent to local 2,3-integrable G'L(2)-structures of 
degree 4, and both objects admit a geometric, coordinate-free classification by the singular foliation 
of shown in the figure. 



Lemma 5.1 seems to be a miraculous coincidence. The bluntness of this relationship between v 



and J{v) prompted me to investigate relationships between the roots of v and the structure of the 



leaf ), yielding this project's main result. Theorem 5.3 It appears that no such relationship 



holds for 2,3-integrable GL(2)-structures of degree n > 5, even though versions of Theorem 4.3 



and Lemma 4/7 exist in those cases. In general, intransitive groupoids and pseudo-groups are very 
poorly understood, and it is generally impossible to explicitly write down the integral manifolds of 
any given singular distribution. If it were not for the coincidence in this case, the leaf-equivalence 
classes would have little utility in understanding the Hessian hydrodynamic PDEs. 

Of course, Hessian hydrodynamic PDEs on m : M'^ — t- M are not the only integrable PDEs of 
interest in mathematics and physics. Two generalizations are important to consider: 

Ql Integrability should be a contact-invariant property of a PDE. The Hessian-only form 
M12, Mi3, M225 ^23, M33) = is not preserved under the full family of contact trans- 
formations. However, the associated GL(2)-structures apparently relied on this form and 



28 



A. D. SMITH 




Figure 2. The leaf-classification of all {B, M,p)2^3- The number of sides on a node 



is the dimension of the bundle after the symmetry reduction from Section 5.1 Rep- 
resentative PDEs for the shaded nodes are included in Section |6l 



its associated CSp{3) transformations. How can this classification be extended to second- 
order PDEs in three independent variables that also include lower derivatives? How does 
the GL{2) geometry generalize to these PDEs? The observations and computations of some 
recent articles may prove very useful [7] [T3] . 
Q2 Very few integrable PDEs are known to exist in more than three independent variables, but 
equations of the form F{uii, . . . ,unn) =0 can sometimes yield distributions of rational nor- 
mal cones of degree n on hypersurfaces M" = F~^{0) C Sym^(M'^) with n = ^A^(A^+1) — 2. 
Results similar to those in Section |4] are known for 2,3-integrability for degrees 5 through 
20, but the foliation by groupoid orbits of V„+4 is not understood, and fc-integrability is 
extremely restrictive [27]. What can integrable GL(2) geometry say about the existence of 
integrable PDEs in more variables? 

On a more detailed level, it would be interesting to study the foliation that appears in the 
present case. Despite the significant computational difficulties, it is important both to produce 
representative PDEs for each root type and to find the root types of the well-known Hessian 
hydrodynamic PDEs. For example, if one can produce Hessian hydrodynamic PDEs as well as 
contact-equivalent PDEs that involve lower-order terms, such computations could open the door 
to Ql, above. Additionally, the lack of surjectivity of T : S — )■ Oj{B) is irritating. What is the 
exact relationship between two PDEs that are leaf-equivalent but do not have overlapping torsion? 
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Finally, the relationship between W, which describes bi-secant surfaces in M, and E/, which 

arises from a planar hyperbolic PDE, is worth pursuing. What is is the nature of the correspon- 
dence? Can every hyperbolic planar PDE appear this way? Can this correspondence provide any 
new information about the Riemann invariants of the Hessian hydrodynamic PDE? 



Appendix A. The Matrix J(T) 

For reference, here is the matrix J{T), listed by column. Note! In the cu columns, the common 
factor of 9216 = 2^°3^ has been removed for clarity. 
The 9216 oj-^ column: 



/ 280T_8T4 - 280T_6T2 \ 

-245T__4T2 + TOT.gTe + 175T_6T4 
70T_4T4 - 210T_2T2 + ISOT.eTg + lOT.gTg 
-I75T0T2 - 35T_2T4 + 35T_qTs + IJbT.^Te 
84r_4T8 + 196T_2T6 - 140T| - l lOToT^i 
-35OT2T4 + 175ToT6 + 175T_2T8 
35OT0T8 - 350T2 
7OOT2T8 - 7OOT6T4 

\ -i40or| + i40or8T4 / 



The 9216 w'^ column: 



/ -1400T_8T2 + 1400T_6To \ 

-385T_8T4 - 840T_6T2 + 1225T_4To 
-280T_4T2 + 1050T_2To - TOT.sU - 700T_6T4 
-910T_4T4 + 280T_2T2 + 875T2 - 240T_6T6 - 5T_8T8 
154070^2 - 952T_2T4 - 28T_6T8 - 560T_4T6 
-105T_4T8 - 1120T_2T6 + 1400T| - I75T0T4 
2IOOT2T4 - I75OT0T6 - 350T_2T8 

2450r| - 1050T0T8 - 1400T2T6 

V -2800^2^8 + 280076^4 / 



The 9216 column: 



/ -2800T_6T_2 + 2800T_8To \ 

-2450T_4T_2 + 875T_8T2 + 1575T_6To 
210T_8T4 - 2100X^2 + 1540T_6T2 + 350T_4To 
1890T_4T2 - 2625r_2To + SbT^sU + 700r_6T4 
1568T_4r4 + 336r_2r2 - 2100r| + 192r_676 + 4718^8 
-2625T0T2 + 1890T_2T4 + 35T_qTs + 700T_^Tq 
210T_4T8 + 1540T_2T6 - 2100T| + 35OT0T4 
-245OT2T4 + I575T0T6 + 875T_2T8 

\ 280oroT8 - 280or2r6 / 
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The 9216 u? column: 
/ 



2800r_4T_6 - 2800T_2T_8 
-1400T_6r_2 - 1050T_8To + 2450X^4 
2100T_4T_2 - 350T_8T2 - irSOT.gTo 
-105T_8T4 + 14007^2 - 1120T_6T2 - M^T.^Tq 
-952T_4T2 + 1540T_2To - 1^T_^T^ - 560T_eT4 
-910T_4T4 + 280T_2T2 + 875T^ - 240T_6T6 - 5T_8T8 
IO5OT0T2 - 280T_2T4 - TOT^gTs - TOOT^^Tq 
I225T0T4 - 385T_4T8 - 840r_2T6 

i4ooror6 - i40or_2r8 



V 



The 9216 co'^ column: 



i4oor_4T_8 - i40or2g 
-700T_4T_6 + 7oor_2r_8 



The A, ip-2, (fio, and ip2 columns: 



-SSOT^^ + 350T_8To 
-350T_4T_2 + 175T_8T2 + ITST.gTo 
84T_8T4 - 140X^2 + 196T_6r2 - 140T_4To 
-35T_4T2 - 175T_2To + 35T_sTe + n^T^&T^ 
70r_4T4 - 210T_2T2 + 130T_6T6 + IQT.^T^ 
70T_6T8 - 245T_2T4 + 175T_4r6 



/T_8\ 
T_4 

To 
T2 



I 



/-16T_6\ 
-14T_4 
-12T_2 
-lOTo 

-8T2 

-6T4 

-4T5 

-2T8 

V J 



/16T_8\ 

i2r_6 

8T_4 
4T_2 


-4T2 

-8T4 
-12Tp 



/ \ 

2r_8 

4T_6 
6T_4 
8T_2 
lOTo 
I2T2 
I4T4 
Vl6Te 
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